Momentum distributions of a-particles from decaying low-lying 12 C-resonances 
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The complex scaled hyperspherical adiabatic expansion method is used to compute momen- 
tum and energy distributions of the three a-particles emerging from the decay of low- lying 12 C- 
resonances. The large distance continuum properties of the wave functions are crucial and must be 
accurately calculated. We discuss separately decays of natural parity states: two + , one 1~, three 
2 + , one 3~, two 4 + , one 6 + , and one of each of unnatural parity, 1 + , 2~, 3 + , 4~. The lowest natural 
parity state of each J n decays predominantly sequentially via the 8 Be ground state whereas other 
states including unnatural parity states predominantly decay directly to the continuum. We present 
Dalitz plots and systematic detailed momentum correlations of the emerging a-particles. 

PACS numbers: PACS: 21.45.-v; 21.60.Gx; 25.70.Ef; 27.20,+n. 



I. INTRODUCTION 

The low- lying resonance states of 12 C have been stud- 
ied over many years both theoretically and experimen- 
tally, motivated partly by their astrophysical importance 
0, 12, S, S H, IE 0, U, i, El- Surprisingly, many issues 
are still not really understood, e.g. the energies, angu- 
lar momenta, structure and decay properties of the res- 
onances. Completely open questions still remain on the 
2 + resonances. Morinaga conjectured in the fifties that 
a 2+ state should exist around 9 MeV as a member of 
the rotational band with the + resonance at 7.65 MeV 
as band-head [ 1| S everal experiments recently provided 
new results [Tllll2l. [l3| but unfortunately no agreement 
has yet been reached for the position and width of the 
first 2 + resonance. 

Attempts to obtain information about the spectrum 
from decay measurements immediately face the problem 
that only the final state is observed. Properties of the 
initial state must then be reconstructed from the momen- 
tum distributions of the three fragments after the decay. 
Both initial state and the intermediate paths connecting 
initial and final states are not observables. These config- 
urations can therefore only be described through model 
interpretations. This is somewhat different in reaction 
experiments, where information can in addition be ex- 
tracted from properties of outgoing particles in transfer 
or scattering reaction. 

If we assume that the initial state is a resonance pop- 
ulated one way or another, and that its decay is inde- 
pendent of the previous history. This is a simplification 
decoupling the formation from the decay in analogy to 
compound nuclear reactions. The decay process can then 
be viewed as a stationary wave function connecting initial 
and final states through a continuous series of intermedi- 
ate configurations. This is equivalent to a time dependent 
process where the initial state, formed at small distances, 
evolves through the intermediate configurations and re- 



sults in the final state at large distances. This implies 
a steady state outgoing flux described precisely by the 
stationary resonance wave function. 



Thus the resonance wave function can be interpreted 
as reflecting the decay mechanisms. Two principally dif- 
ferent modes are traditionally considered, i.e. sequential 
decay via an intermediate two-body configuration, and 
decay directly into the three-body continuum. In both 
cases the final state is embedded in the three-body con- 
tinuum and the modes can only be distinguished if the 
momentum distributions carry unique information char- 
acterizing one of the modes. Otherwise the distinction 
becomes fluent or a matter of an artificial, although per- 
haps more precise, model definition. Previous approaches 
to describe this type of observables have been performed 
mainly for the 1 + states, e.g. Faddeev calculations 
R-matrix computations, which describe their decay as 
sequential 15], and Kurchatov fitting, which describes 
them as direct or democratic [16j . 



The purpose of this paper is to present a-particle 
momentum distributions and Dalitz plots 1171 after de- 
cays of all the computed 12 C-resonances [18| below the 
proton separation threshold at an excitation energy of 
15.96 MeV, where only 3a-decay is possible. These dis- 
tributions should help to establish spins and parities of 
the yet unknown levels. They provide then information 
about structures of initial and intermediate states. Com- 
bined with the measurements a more complete picture 
of the 12 C-spectrum and the decay mechanisms should 
then emerge. In section 2 we first sketch the theoretical 
framework and the choice of interactions. The results are 
presented and discussed in section 3 for both unnatural 
and natural parity states. Section 4 contains a summary 
and the conclusions. 
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II. THEORETICAL FRAMEWORK 

The resonances decay into three particles, therefore we 
need a theoretical tool to describe this three-body con- 
tinuum structure. We employ the hypersphcrical com- 
plex rotated [H, H3| adiabatic expansion [2l| in coor- 
dinate space to compute bound states and resonances. 
This method is able to deal with several simultaneously 
bound and nearly bound two-body states in different sub- 
systems. Relatively large distances can often be calcu- 
lated accurately with a specific choice of basis and par- 
tial waves. The Fourier transform of the wave function 
provides the observable momentum distributions. The 
three-body model consisting of a-particles requires inter- 
actions which reproduce energies and scattering proper- 
ties of the a-a system. 

A. Practical procedure 

We describe 12 C as a 3a-cluster system at all dis- 
tances. We use Faddeev equations and solve them in 
coordinate space using the adiabatic hyperspherical ex- 
pansion method [2(| HJ, [22j • The hyperspherical coordi- 
nates consist of the hyperradius p and five generalized 
angles. The angular Faddeev decomposed wave func- 
tions, $„ JA f = Ei=l $ n 

j M , are chosen for each p as the 
eigenfunctions of the angular part of the complex scaled 
(r — ► fcxp(i9)) Faddeev equations 

(T n - X n ) + 5pV ! $ n ,/Af=0 i = 1, 2, 3 . (1) 
n 

Th is the angular part of the kinetic energy operator 
and Vi is the potential between particles j and k, be- 
ing {i,j, k} a cyclic permutation of {1,2,3}. The total 
wave function, ^ JM } is expanded on the hyper-angular 
eigenfunctions, i.e. 

* JM = 4iE/»0 9 ) <i> »^ ) ' ( 2 ) 

P n 

where the p-dependent expansion coefficients, f n (p), are 
the hyperradial wave functions obtained from the coupled 
set of hyperradial equations 

/ d 2 15/4 2m r . , . , A 

( -Qf + ~jr + -ffl W n{p) + VM - E]j f n (p) 

OG 

= Y, P nn'fn'(p) (3.) 
n' = l 

W n (p) are the angular eigenvalues of the three-body sys- 
tem Hamiltonian with fixed p, V$b is the three-body po- 
tential, E is the three-body energy and P nn > are the non 
adiabatic terms. The eigenvalues W n (p) of the angular 
equations eq. JT]) serve as effective potentials. 

In order to obtain the resonances we use the complex 
scaling method. According to this method, the energy 



and width of a resonance state are associated with the 
complex eigenvalues of a certain analytically continued 
Hamiltonian operator. The appropriate operator results 
from the rotation of the position vectors of the ordinary 
Hamiltonian into the complex coordinate plane 

f->re l9 0>O,real. (4) 

This gives rise to the complex-rotated Hamiltonian 

H e (r) = H(re ie ) . (5) 

The complex energy of a resonance corresponds to 
a pole in the momentum-space wave function, while 
in coordinate-space this form corresponds to a large- 
distance asymptotic wave function consisting of outgoing 
waves. In other words, the three-body resonance corre- 
sponds to a complex energy solution Eq = Er — i Ej of 
the system ^ with the asymptotic boundary condition 
of an outgoing wave in every channel n 

f n {p^^) = C n e +l ^ , (6) 

where C n is an asymptotic normalization coefficient and 

k = \J 2mE/h 2 is the three-body momentum or the con- 
jugate of p. It has been seen that this boundary condition 
determines that the scattering matrix has a pole at the 
complex energy Eq, being Er the position of the reso- 
nance and r = 2Ei its width. 

The 12 C-resonances are not necessarily of three-body 
character even though this by definition must be the case 
at large distances for 3a-decay. We use the three-body 
model also at small distances because, like in Gamow's 
theory of a-decay, the detailed structure at small dis- 
tances is not important for the decay properties which 
only require the proper description of the emerging three 
particles. We use the three-body short-range potential to 
adjust the corresponding small-distance part of the effec- 
tive potential to reproduce the correct resonance energies 
which are all-decisive for decay properties as evident in 
the probability for tunneling through a barrier. 

At intermediate distances the three a-particles are 
formed and the potential has a barrier that determines 
the partial width of the resonance. At large distances 
the resonance wave functions contain information about 
distributions of relative energies between the three par- 
ticles after the decay. These properties are connected 
to the many-body properties at small distances via pre- 
formation factors, as in a-decay. An adjustment of the 
resonance energy is then needed. After complex rota- 
tion the resonance wave function is characterized by an 
exponential fall-off at large distance. Thus the crucial 
information is found in relative sizes of the very small 
values, /„, of the resonance at large distances which are 
very difficult to compute accurately especially when the 
Coulomb interaction is present. 
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B. Momentum distributions 

The complex scaled coordinate space resonance wave 
function should be rotated back to real coordinates and 
Fourier transformed to provide the observable momen- 
tum distributions. Unfortunately the corresponding in- 
tegral is not convergent and a regularization procedure 
has to be applied. The origin is simply that the result- 
ing wave function should be a non-normalizable outgoing 
plane wave at large distances. We overcome this prob- 
lem with the Zeldovic regularization procedure which is 
well defined for short-range interactions [20]. In total 
this amounts to using the angular part of the coordinate 
space wave function at a large hyperradius, but inter- 
preted as the momentum space wave function. Inclusion 
of the Coulomb interaction is achieved by treating it as an 
ordinary potential up to a large value of the hyperradius 
and then extrapolate the diagonal parts of the adiabatic 
wave functions with the numerically obtained Coulomb 
and centrifugal potentials. 

Two different cases must be treated, i.e. sequential and 
direct decays distinguished theoretically by the structure 
of the adiabatic wave functions [23| . Direct decay is char- 
acterized by structures where all particles are far apart 
and as the hyperradius increases all distances increase 
proportionally. The Zeldovic regularized Fourier trans- 
form of the resonance wave function gives in this case 
the momentum distributions [2~jj . 

Sequential decay is characterized by a wave function 
describing a bound-state like structure of two close-lying 
particles supplemented by the third particle far away. For 
a complex scaled wave function such a structure would be 
that of a two-body resonance, provided the rotation angle 
9 is larger than the angle corresponding to the energy 
and width of this resonance. These structures approach 
two-body bound state configurations as the hyperradius 
increases. 

However, Fourier transformed and rotated back to the 
real axis, the wave function should at large distances 
approach the description of the third particle (plane or 
Coulomb wave) leaving the decaying resonance which has 
the given two-body energy and width. This is two se- 
quential two-body decays, hence the characterizing nota- 
tion. The resulting momentum distributions cannot be 
obtained from the rotated wave function but should in- 
stead be calculated from the correct physical description 
of two two-body decays. This results in a Breit-Wigner 
distribution for the third particle with a width equal to 
the sum of two-body and initial three-body resonance 
widths peaking around the energy found by subtracting 
the two-body from the three-body energies. 

C. Interactions 

The basic ingredients are the two-body interactions,!^, 
between particles j and k, where k} is a cyclic per- 
mutation of {1,2,3}. First Vi must reproduce the low- 



energy two-body scattering properties which can be ob- 
tained independently for each partial wave resulting in 
angular momentum dependent or non-local interactions. 
We rely on the experience gained previously especially 
through [l8| , and we choose an Ali-Bodmer potential [25] 
slightly modified in order to reproduce the s-wave reso- 
nance of 8 Be. The phase shifts are essentially unchanged 
and reproduce a — a scattering data but in order to de- 
scribe sequential decays properly the two-body subsys- 
tems must also have the correct energies. In total we use 
a potential given as 

V aa = (l25P i=0 + 20P ;=2 ) e - r2 /i.53 2 

-30.18e- r2/2 ' 852 , (7) 

where lengths are in fm and strengths are in MeV. The 
operators Pi project on angular momentum. 

The three-body resonance energy and wave function 
can now be computed but the energy usually does not 
coincide with the measured value. It may be close, indi- 
cating that the three-body structure is nearly correct. 
Then only fine-tuning is needed due to the neglected 
smaller three-body effects of polarization or excitations 
of intrinsic particle degrees of freedom or off-shell effects. 
We emphasize that only three-body effects are missing 
since the two-body data already is reproduced by the 
phenomenological two-body interactions. We then cor- 
rect the energy by including a diagonal three-body short- 
range interaction chosen to be Gaussian in hyperradius, 
i.e. Vsb = Sexp(— p 2 /b 2 ). The structures of the reso- 
nances are then maintained [26]. A larger range corre- 
sponding to a third order power law is not selected as e.g. 
in [2jJ where it is used to compensate for the limitation 
in Hilbert space due to the hyperharmonic expansion in 
only one Jacobi coordinate. Our better basis confines 
the three-body interaction to be genuinely of short-range 
character. 

In the actual parameter choice we prefer to maintain 
the same values of b and S for different states with the 
same angular momentum and parity J 71 but allow vari- 
ation with J 71 '. To see the systematic behavior we then 
decided to fix b = 6 fm corresponding to the hyperra- 
dius obtained when the three alphas are touching in an 
equilateral triangle. The strength S is then adjusted to 
reproduce one of the observed resonance energies. The 
main dependence is indirect through the variation of the 
three-body energy and much less through the shape of the 
total potential 18]. In this way we attempt to separate 
the effects of the initial many-body structure from the 
symmetries related to the angular momentum conserva- 
tion. The strongest influence is expected from Coulomb 
potentials and centrifugal barriers. 

III. COMPUTED DISTRIBUTIONS 

We find 12 C-resonances below 15.96 MeV for most an- 
gular momenta J < 6 and all parities, i.e. two + , three 
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2+, two 4+, and one of each of 1 ± , 2~, 3*, 4~ and 6+ 0. 
Their structures were described in detail in [18] including 
the variation with possible interaction parameters. How- 
ever, only small and intermediate distance properties are 
important for energies, widths and partial wave decom- 
position. The final state momentum distributions after 
decay arise from the large distance properties which are 
much more difficult to determine numerically. 

The procedure is to compute ratios of radial wave func- 
tions at large distances. This supplies the relative weights 
on the contributions from each of the adiabatic wave 
functions. First we have to remove the contributions 
from the wave functions corresponding to population of 
two-body resonances. These fractions must be computed 
as consecutive two-body decays and their contributions 
added to the remaining results from direct decays which 
are found by absolute square of the wave function at a 
large hyperradius followed by integration over the unob- 
served angular variables. 

The asymptotic large-distance behavior should be 
reached by increasing the partial waves and the basis 
size used. This convergence can be tested by showing 
independence of the results with variation of the largest 
value of the hyperradius. Failing the test implies that the 
basis size is too small, or contrarily a larger hyperradius 
can be compensated by a larger basis producing the same 
result. It is then economical to get stability for a hyper- 
radius and basis as small as possible. In most cases we 
find that the asymptotic behavior is reached for hyper- 
radii larger than about 60 fm. There is a small variation 
of the distributions from 70 to 100 fm, and we have cho- 
sen 80 fm as the value of p where the energy distributions 
are computed. 

The results fall in two groups of natural and unnatural 
parity states, e.g. implying that sequential decay through 
8 Be(0 + ) is either allowed or forbidden by conservation of 
angular momentum and parity. Decay through 8 Be(2 + ) 
is possible in both cases but this state is rather broad and 
the result would be hard to distinguish from direct decay. 
We see no indication of population of this channel in 
the numerical results. To optimize the accuracy we then 
maintain as small a scaling angle as possible consistent 
with distinct separation of the three-body resonance from 
the background continuum. 



A. Unnatural parity states 

These states are 1 + , 2 _ , 4~ and 3 + and our basis de- 
scribes them as decaying directly although analyses of 
measured distributions employ interpretations as sequen- 
tial through 8 Be(2+) [11 [H 111 • 

The lowest 1 + state was briefly discussed previously in 
[H,[3(|. Experimentally two 1 + states, isospin and 1, 
are known but we find only one reflecting that we are con- 
fined to isospin by using a-particles as building blocks. 
Both states are very far from resembling a-cluster states. 
Still the decays of both states must proceed through the 



same a-cluster configurations, although the weights on 
the adiabatic potentials might differ from state to state. 
Underlying many-body effects are beyond the present 
model but we can pinpoint the neglected effects, i.e. the 
preformation factors established at small distances where 
the many-body problem is constraint into a three-body 
problem, and better three-body potential to account for 
the transition between the N- and three-body degrees of 
freedom at short and large distances, respectively. For 
these reasons the contributions from the individual adi- 
abatic potentials could differ for decays of these two 1 + 
states of different isospin. 

We first focus on the isospin 1 state at an excitation 
energy of 14.98 MeV (7.70 MeV above threshold). We 
adjust the three-body potential and compute the energy 
distributions shown in Fig[TJ The upper part exhibits 
the Dalitz plot and the lower part projects the distribu- 
tion on the axis with one a-particle energy. The latter 
is computed by using Monte Carlo integration over all 
phase space directly from the wave function. 

The measured distributions [3l| are very uncertain first 
of all because the lower-lying isospin zero 1 + state at 
12.70 MeV (5.42 MeV above threshold) also is populated 
via feeding from a gamma transition between the two 
1 + states. This contribution is not easily removed from 
the existing data to allow a clean comparison. A better 
analysis or a new experiment measuring the a-decay of 
the T=l 1 + state in complete kinematics is required. Our 
computed result is almost identical to the distributions, 
measured and calculated, for the 12.70 MeV state [23ll29j 
if the difference in available energy is corrected for. The 
distributions in FigfT] are then direct prediction based on 
the assumptions that the isospin zero components in both 
states are equally populated and decay through the same 
mechanism. 

A test of this prediction would provide interesting in- 
formation about the dynamics of isospin mixing. Two 
extremes can be imagined, i.e. the same isospin com- 
ponents can be present from small to large distance re- 
sulting in the same distribution, or different complicated 
many-body structures at small distances clusterize into 
a-particles around the nuclear surface and proceed to 
detection at large distances. We know that the partial 
decay widths for both states are much smaller than pre- 
dicted from the cluster model [H| but this information 
does not prohibit the momentum distributions from be- 
ing almost identical. 

In the computation we find only one of each state of 
even J and negative parity, i.e. one 2~ one 4 _ state. The 
experimental spectrum has two states of 2~ where the 
highest at 13.26 MeV (5.98 MeV above threshold) only 
is tentatively assigned to have 2~ [32j while no 4~ state is 
found experimentally. It is then tempting to believe that 
this state at 13.26 MeV really is a 4~ state as indicated by 
our computations [33j. This new sp in-p arity assignment 
has also been suggested recently in [34j |. 

One way to decide which spin and parity is correct is 
to measure the momentum distributions of the fragments 
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FIG. 1: Dalitz plot (upper part) and the a-particle energy distribution (lower part) for the (l + ,2 ,3 + ,4 )-resonances at an 
excitation energy of (14.98, 11.80, 14.40, 13.26) MeV or (7.70, 4.53, 7.13, 5.98) MeV above the 3a threshold, which is 7.275 MeV 
above the ground state. We have performed a Monte Carlo integration over the phase space, which due to the statistical nature 
produces the unphysical fluctuations. 




FIG. 2: The angular distributions of the directions between two particles and their center of mass and the third particle for 
the (l+,2-,3+,4-) -resonances in Fig[T] We have performed a Monte Carlo integration over the phase space, which due to the 
statistical nature produces the unphysical fluctuations. 



emerging after decay. Usually this carries distinct signa- 
tures of the angular momentum of the decaying state. In 
[35l [3fjj it is shown that, even within non-sophisticated 
theoretical models, the basic signatures of the angular 
momentum are present in the experimental data. We 
first turn to the energy distributions in Fig[T] Both 2~ 
and 4 _ show very similar distributions but the peaks ap- 
pear at slightly higher values for the 4~ state. However, 
the two-dimensional Dalitz plots differ more from each 



other. Both have the triangular symmetry but the 2~ 
resonance have virtually nothing in between these peaks 
in contrast to the somewhat more smeared out distribu- 
tions of the 4~ resonance. 

In the computations we find also a 3 + resonance for 
which there is no experimental evidence, but it has been 
suggested in [34| to assign these quantum numbers to 
the state at about an excitation energy of 13.35 MeV. 
Theoretically a 3 + state has also been found in 0] ■ With 
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a reasonable three-body strength, -20 MeV, placing the 
state at 14.40 MeV (7.13 MeV above threshold), we find 
the energy distributions in fig. [T] The distribution is very 
broad but peaked at intermediate energies. This is seen 
to arise from a Dalitz plot distribution with a small hole 
in the middle surrounded by a close-lying dense circle 
and a much larger diffuse distribution. 

The angular momentum may leave an even more dis- 
tinct signature in the angular distributions, shown in 
Fig|2j of the directions between two particles and their 
center of mass and the third particle. The information is 
then directly about the corresponding angular momen- 
tum denoted e. the angular momentum of the 
third particle relative to the center of mass of the other 
two with the relative angular momentum t x . We see that 
the angular distribution patterns are quite different for 
different states. The 1+ distribution shows two broad 
peaks separated by a minimum with vanishing probabil- 
ity at an angle of tt/2. 

This reflects that the partial wave components in the 
angular wave function are a linear combination of only 
(4, 41 = (2, 2), (4, 4) each coupled to the resulting value 
of 1 [23J. Choosing the specific directional angles of 
4>x = <t>y — @x — and 9 y = tt/2 wc find that only pro- 
jection quantum numbers of m x — and m y = 0, ±2, ±4 
give non- vanishing contributions. This is only consis- 
tent with a projection of the total angular momentum 
M = m x + rriy since M = m y = ±1 gives zero. How- 
ever when all projections are zero the Clebsch-Gordan 
coupling coefficient is also zero. The observables in FigJ^] 
reveal information about the intrinsic angular momenta 
used to construct the wave function. 

In contrast both 2~, 3 + and 4~ have peaks in the dis- 
tributions at 7r/2. The different shapes can be traced 
back to the different partial wave decomposition com- 
puted and discussed in [23], i.e. 2~ has about 40% to 
60% of l y — 1, 3, while 4~ is dominated by £ y = 3, and 
3 + has about twice as much £ y — 2 as £ y — 4. These fea- 
tures are clearly distinguishable demonstrating that these 
observables can be used to determine the large-distance 
structure of these resonances. The initial state can still 



only be determined through the theoretical information 
about the dynamical evolution of the resonances. 

The one-dimensional distributions in Fig[T]can be used 
to extract the distributions of how far the three particles 
are from each other 37]. This is visualized by a triangle 
with a particle in each corner moving apart from their 
common center of mass. In particular the distributions 
of the ratio, x/y, of the distances between two particles, 
x, and their center of mass and the third particle, y, are 
shown in Fig|3J Since all Jacobi systems are identical we 
do not have to distinguish between Jacobi sets. Unfor- 
tunately the symmetric wave function then do not allow 
distinction between these identical particles. 

With several peaks as for the 1 + resonance the inter- 
pretation is obvious, namely that each peak contains one 
a-particle. The triangular geometric structure for the 
decay of this isospin one 1 + state then corresponds to 
side ratios of 2.2:1.8:1 of the triangle. For the other un- 
natural parity states only one broad peak is seen close to 
the value 1. For an equilateral triangle the x/y-ratio is 
2/ y/3 ss 1.15 which then is the only value where a narrow 
peak is possible. Otherwise a broader peak must cover 
overlapping distributions deviating somewhat from the 
equilateral triangle and corresponding to similar but less 
symmetric configurations. 



B. Natural parity states 

These states are 0+, 1", 2+, 3", 4+ and 6+. They 
can decay via the energetically favorable 8 Be(0 + ) which 
asymptotically must be described by one of the adiabatic 
potentials with the 8 Be+a structure. The signature is 
simply that this potential approaches the complex en- 
ergy of the 8 Be(0 + ) resonance. The radial resonance 
wave functions at large distances determine the popu- 
lation fractions for each of the adiabatic potentials. In 
particular we can find the fraction of decay proceeding 
sequentially through this + -state, and furthermore we 
can compute the related distributions as two consecutive 
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FIG. 4: Dalitz plot (upper part) and the a-particle energy distribution (lower part) for the (Oj, 2^~, 2jj~)-resonances at an 
excitation energy of (11.22, 11.76, 13.76) MeV or (3.95, 4.48, 6.49) MeV above the 3a threshold, which is 7.275 MeV above the 
ground state. We have performed a Monte Carlo integration over the phase space. The sequential decay, respectively of 59%, 
15%, 4%, through 8 Be(0+) is removed. We label as in table U 



12 C(1 ) 12 C(4 + 2 ) 12 C(6 + ) 
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FIG. 5: Dalitz plot (upper part) and the a-particle energy distribution (lower part) for the (1 , 4^~, 6 + )-resonances at an 
excitation energy of (10.88, 14.10, 14.40) MeV or (3.61, 6.83, 7.13) MeV above the 3a threshold, which is 7.275 MeV above the 
ground state. We have performed a Monte Carlo integration over the phase space. The sequential decay, respectively of 70%, 
20%, 5%, through 8 Be(0 + ) is removed. We label as in table U 
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TABLE I: Energy above the triple-Q threshold, excitation en- 
ergy and estimated amount of sequential via 8 Be(0 + ) and di- 
rect decays for the natural parity states of 12 C. If necessary we 
label the resonances with increasing energy above threshold. 



J 71 " E QQC1 (MeV) E exc (MeV) sequential direct 



3- 

4 
1- 

0+ 

2+ 
6+ 



0.38 


7.66 


95% 


5% 


1.38 


8.66 


97% 


3% 


2.33 


9.60 


96% 


4% 


3.25 


10.52 


92% 


8% 


3.61 


10.88 


70% 


30% 


3.95 


11.22 


59% 


41% 


4.48 


11.76 


15% 


85% 


6.49 


13.76 


4% 


96% 


6.83 


14.10 


20% 


80% 


7.13 


14.40 


5% 


95% 



two-body decays. 

The result is one peak close to an energy of E max = 
2E a /3 with a width roughly equal to the width of the de- 
caying state, and a broader square-like peak at an energy 
of about E max /A determined by kinematics. Here we as- 
sumed vanishing energy and width of 8 Be(0 + ), otherwise 
the peak positions and widths should be modified. The 
Dalitz plots should also reflect these features by show- 
ing one high-energy, almost vertical, single-a distribu- 
tion, and two separated (for each of the other a-particles) 
more horizontal distributions corresponding to a broader 
peak after projection on the single a-energy x-axis. 

The angular distribution from sequential decay 
through 8 Be(0 + ) must reflect the behavior of the angu- 
lar momentum £ y precisely as for ordinary decays of a 
quantum state of given angular momentum. The direct 
decay is expected to give a relatively broad distribution 
shifted from the central value at half the maximum en- 
ergy by an appropriate average over the combinations of 
angular momentum phase space factors. This can also 
be interpreted geometrically as an expanding triangular 
configuration with given side ratios. 

First we extract the percentage of sequential decay via 
8 Be(0 + ) and direct decay for the natural parity states, 
see table |H The lowest-lying natural parity states of each 
J n (0 + , 2+, 3~ and 4+ states with excitation energies 
7.66 MeV, 8.66 MeV, 9.60 MeV and 10.52 MeV) seem 
to be completely dominated by decays via 8 Be(0 + ). In 
contrast, the highest-lying 2 + state at 13.76 MeV excita- 
tion energy and the 6 + state at 14.40 MeV excitation en- 
ergy only have small fractions decaying through the 8 Be 
ground state. In the remaining cases (1 — , + , 2 + and 4 + 
states with excitation energies 10.88 MeV, 11.22 MeV, 
11.76 MeV and 14.10 MeV) both mechanisms are com- 
parable. 

The lowest of the two + resonances is the so-called 
Hoyle state, which plays an important role in nuclear 
astrophysics. According to our computation, it decays 
almost entirely sequentially. Very little is left for the di- 
rect decay which therefore is not shown. The experimen- 
tal distribution is also consistent with complete domina- 



tion of sequential decay as in our computation (23[. We 
also omit the other three natural parity resonances dom- 
inated by sequential decays. We concentrate instead on 
the 6 resonances where a substantial amount is direct de- 
cay. These distributions are shown in Figs. [4] and \5\ after 
removal of the contributions from the sequential decay 
through the 8 Be ground state. The experimental anal- 
yses can extract the trivial contribution from the decay 
through the 8 Be ground state. It is therefore straightfor- 
ward to make a comparison with the experiment. Both 
Dalitz plots and projected single-a energy distributions 
are shown. 

The higher-lying + resonance has a large width of 
about 3.5 MeV. On top of this difficulty the population 
through beta-decay of the corresponding energy region 
leads to violation of the independent approximation of 
formation and decay of the resonance. The main effect 
is a shift in energy of the resonance position. In any 
case this state has a significant probability of decaying 
directly into the three body continuum. This part, shown 
in fig. HI exhibits a triangular structure in the Dalitz plot, 
but now we find one low-energy cv-particle and two of 
moderate energies. This is in almost complete contrast 
to the sequential decay where one energy is high and two 
are small. 

Next we focus on our results in connection with the 
existence and position of low-lying 2 + resonances which 
still is an open question for the 12 C nucleus. The old sug- 
gestion is that the Hoyle state should be the band-head 
followed by a 2 + state at around 10 MeV [l[. There are 
experimental indications for the existence of such a state 
[29| | but no consensus has so far been reached. On the 
other hand other theoretical models, also cluster models, 
find three 2 + resonances in this energy region [38[ ; in 4 1 
two 2 + excited states are found in this region, while in 6| 
one 2 + state appears below 12.3 MeV excitation energy. 
We find rather different structures for these three states, 
still all of cv-cluster structure. Each of them is dominated 
by its own adiabatic wave function corresponding to three 
different low-lying adiabatic potentials with differing par- 
tial wave decomposition [l8| • Most likely these states are 
hidden behind broad states of roughly the same energy. 
They are therefore extremely difficult to distinguish from 
the background in any of the experiments. 

Their decay properties also vary substantially, e.g. the 
percentage of sequential decay through 8 Be(0 + ), see table 
[H The lowest state almost exclusively decays sequentially 
while the other two mostly decay directly. In fig. [4] we 
see that the direct parts give very broad distributions. 
For the second 2 + resonance all three cv-particles emerge 
with large probability with similar kinetic energies. For 
the third 2 + resonance the distribution is more diffuse 
and the energies are more unevenly divided resulting in 
a structured but relatively broad distribution. 

We continue with the 1~ state, both Dalitz plot and 
one-dimensional projection are shown in Fig. [5] for the 
30% decaying directly into the three-body continuum. A 
similar triangular structure as for the second + reso- 



9 




Jt/2 71 Jt/2 t 



FIG. 6: The angular distributions of the directions between 
two particles and their center of mass and the third parti- 
cle for the (0^, 1~, 2^, 2^~, 4j, 6 + )-resonances in Figs(4]and[5] 
We have performed a Monte Carlo integration over the phase 
space. The sequential part is removed as in figs. |4]and[5] We 
label as in table [I] 



nance is seen although substantially more smeared out 
resulting in two overlapping broad peaks after projection 
on the x-axis. 

Both the 3~ resonance and the lowest of the two 4 + res- 
onances are almost completely dominated by sequential 
decay. The second of the 4 + resonances gives a rather 
diffuse distribution of kinetic energy of the ev-particles, 
see fig. [3] It resembles somewhat the distribution from 
the third 2+ resonance except that the small probability 
holes in the Dalitz plot now also are smeared out. This 
distribution is again almost the opposite of the sequen- 
tial decay distribution with one high and two low energy 
particles. The 6 + resonance has a symmetric distribution 
extending about 1 MeV around a central region where all 
energies are roughly equal. 

We now turn to the other type of information found 
in the angular distributions which exhibit the correlated 
directions of emergence. Obviously the sequential de- 
cay through the 8 Be(0 + ) state must be with the third 
a-particle in the opposite directions of 8 Be. The only 
information here is then about the partial wave compo- 
nent (£ y ) of that third particle relative to 8 Be. Angular 




x/y x/y 



FIG. 7: The distributions of the ratio of the distances between 
two particles and their center of mass and the third particle 
for the (0j,l~,2+,2+,4+,6 + )-resonances in FigH The se- 
quential part is removed as in figs. [4] and [5] We label as in 
table |U 



momentum conservation then requires the total angular 
momentum J = £ y . Thus the most interesting new infor- 
mation is contained in the directly decaying parts shown 
in FigEl These distributions also vary from state to state 
reflecting the structure in terms of partial waves as dis- 
cussed in [TH - 

The distribution corresponding to 0^-state is essen- 
tially from the isotropic distribution of £ y — modified 
by a smaller contribution from £ y = 2 with maxima at 
7r/4 and 37r/4 separated by zero probability at 7r/2. The 
distribution corresponding to 1~ shows two peaks sepa- 
rated by a small minimum at tt/2. The largest partial 
waves are here £ y = 1,3. The angular distributions of 
both the second and third 2 + resonance seem to contain 
a narrow peak on top of a broader one. These struc- 
tures are due to large contributions from £ y = supple- 
mented by contributions from £ y = 2 and £ y — 4, respec- 
tively. Finally, the distributions form 4^ and 6 + both 
exhibit one smooth, and for 6 + also relatively narrow, 
peak around tt/2. The partial wave structures of these 
states are mainly £ y = 2, and £ y = 2,4, respectively 

We again attempt to extract the geometric structure 
of the dominating triangular decay configurations. The 
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results for the ratio between one pair of particles and 
their center of mass and the third particle are shown in 
FigEl They are all rather similar with a relatively broad 
peak around 1, but for 1" and 2^ with more structure at 
a larger ratio suggesting another peak. As in Fig[3] the 
peaks must cover overlapping distributions to correspond 
to an almost equilateral triangle. In the case of + a very 
broad peak appears around 3, and the other two peaks 
are around 0.8. This gives rise to an obtuse triangle with 
side ratios 1.7:1:1. 



IV. SUMMARY AND CONCLUSIONS 

We have computed the a-particle momentum distri- 
butions of 14 three-body decaying low- lying 12 C many- 
body resonances with 10 different angular momenta and 
parities. The results are exhibited as single a energy 
distributions and energy correlations of Dalitz plots. We 
assume that the decays of the resonances are independent 
of their formation as for compound nuclear reactions. We 
use a three-a cluster model to describe all states even at 
small distances where the cluster model sometimes fails 
badly and the many-body structure is indispensable for 
a structure computation. The idea is, the same as for 
the classical a-emission, that three a-particles must be 
formed at small or intermediate distances as they emerge 
at large distances after the decay. Thus the small dis- 
tance properties should only supply boundary conditions 
and impose energy and angular momentum conservation. 
This we mock up in the 3a cluster model by a three-body 
interaction adjusted to reproduce the resonance energy. 
Again a simple analogy is found in the preformation fac- 
tors in a-emission. 

An extreme example is the isospin 1 state which can- 
not be formed by a-clusters. Its a-decay width is conse- 
quently very small but still the resulting distributions are 
with the present assumptions predicted to be essentially 
the same as the 1 + isospin state. 

For three-body decays the interest, and the compli- 
cation, is how the energy is shared between the three 
particles. This is determined by the "dynamic evolu- 
tion" of the resonances, i.e. by the change in structure 
from small to large distances. To a large extent the de- 



cisive properties are symmetries from angular momen- 
tum and parity conservation. The resulting momentum 
distributions carry information about both initial reso- 
nance state and the intermediate configurations (decay 
mechanisms). The only energetically allowed two-body 
structure is the ground state of 8 Be. Sequential decay 
through this state is dominating for natural parity states 
for the lowest resonance of a given angular momentum. 
The momentum distributions for the fractions decaying 
directly are predicted for all resonances below the proton 
separation threshold. 

Whenever possible we give a geometric description of 
the parts decaying directly to the three-body continuum. 
This is expressed as side ratios of the a-particles emerging 
in a triangle. 

The Dalitz plots and a-energy distributions differ from 
state to state. A complementary observable is the corre- 
lation between the direction of one particle and the center 
of mass of the other two. These distributions could be 
used to assign spin and parity to these decaying states as 
soon as sufficiently accurate experimental data become 
available. The directly measured angular distribution 
must contain information about the angular momentum 
of one particle with respect to the center of mass of the 
other two particles at large distances. Since several par- 
tial waves may contribute this information is not unique, 
and may have to be supplemented with other informa- 
tion. Furthermore, the uncertainty remains of how the 
measured large-distance properties reflect the small and 
intermediate-distance structures of the resonance wave 
function. Only a theoretical model can provide this con- 
nection. 

In conclusion, we provide systematic and detailed de- 
cay information (fraction of sequential decay, Dalitz 
plots, single-a energy distributions, momentum direction 
correlations), which can be compared to upcoming ex- 
perimental data, for each of the 14 lowest 12 C resonances 
decaying by 3a-emission. 
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